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摘要 : 本 文 基于 Python 设计 和 实现 自动 测 光 程序 。 首 先 扫 描 文件 名 实现 文件 分 


类 ; 然后 进行 数据 检查 、 预 处 理 与 孔径 测 光 ， 并 保存 测 光 数据 ; 最终 由 观测 者 指 


定 目标 星 与 比较 星 , 程序 从 测 光 结果 中 提取 相应 星 等 与 观测 时 间 等 数据 ,获得 基 


于 较 差 测 光 法 的 光 变 曲线 。 为 解决 数据 提取 中 星 图 匹配 问题 , 提出 自动 奶 踪 移动 


目标 的 基于 向 量 特 征 星 图 匹配 算法 , 能 自动 匹配 偏 移 较 大 的 星 图 ， 并 跟踪 不 断 移 


动 的 小 行星 。 使 用 本 程序 分 别处 理 双 星 和 小 行星 观测 数据 ， 并 与 IRAF (Image 


Reducation and Analysis Facility〉 测 光 结 果 比 较 ， 测 光 结 果 最 大 偏差 分 别 为 


0.04mag 和 0.043mag， 标 准 偏差 分 别 为 土 0.005mag 和 土 0.007mag。 
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天 文 观 测 是 获取 天 体 电磁 辐射 流量 的 常用 方法 之 一 。 天 体 的 许多 特殊 现象 如 


双星 撼 食 、 人 恒星 脉 动 、 色 球 活动 和 糙 发 、 超 新 星 爆发 、 伽 马 射 线 暴 等 ， 都 可 以 从 
测 光 中 获得 数据 。 随 着 现代 技术 的 发 展 ， 人 们 普遍 使 用 二 维 探测 器 CCD 测 光 ， 


它 具 有 量子 效率 高 、 波 长 响应 范围 宽 等 优点 。 由 于 绝对 定 标 比 较 困 难 ， 人 们 常用 


较 差 测 光 研究 天 体 光 度 的 变化 ， 以 消除 来 自 热 发 映 和 菊 光 发 射 的 天 空 背景 辐射 。 


近年 来 ， 随 着 地 面 和 空间 巡天 项 目的 开展 (如 国内 的 GWAC 大 视 场 巡 天 站 和 即 
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将 发 射 的 中 国 空间 站 巡 天 空间 望远镜 站 等 )， 海 量 的 测 光 数据 需要 及 时 处 理 。 此 


外 ， 观 测 中 发 现 一 些 瞬 变 源 如 恒星 耀 发 ”和 伽 马 射线 暴 ""， 需 要 适时 进行 数据 


处 理 。 因 此 开发 一 套 自动 测 光 分 析 软 件 十 分 必要 。 


IRAF 是 美国 国家 天 文 台 开 发 的 天 文 图 像 处 理 和 分 析 软 件 ， 广 泛 用 于 处 理 测 


光 数据 " “。 但 该 软件 使 用 繁琐 ， 需 要 较 多 的 人 工 操作 。 为 提高 测 光 效率 ， 有 研 


究 者 以 IRAF 为 基础 开发 自动 测 光 程序 包 ![ ， 在 部 分 望远镜 中 得 到 应 用 ， 但 未 


能 普及 。 随 着 Python 语 


言 的 兴起 ， 出 现 了 处 理 天 文 数 据 的 程序 包 ， 使 得 我 们 能 


够 基于 Python 语言 ， 开 发 出 一 套 自 动 测 光 程 序 ， 可 对 测 光 数据 分 析 并 提取 光 变 


曲线 。 


1. 算法 设计 


为 消除 大 气 满 流 等 天 气 的 影响 , 通常 采用 较 差 测 光 法 获得 目标 光 变 曲线 。 即 


以 同 视 场 中 另 一 颗 亮 度 不 变 的 恒星 为 比较 星 , 计算 目标 星 与 比较 星 的 亮度 差 ， 从 


而 获得 准确 的 光 变 曲线 。 本文 程序 测 光 流 程 分 为 两 部 分 , 对 图 像 中 所 有 恒星 测 光 ， 


根据 需要 提取 目标 星 亮度 。 


1.1 自动 孔径 测 光 
1.1.1 数据 预 处 理 


观测 中 , CCD 噪声 、 本 底 以 及 像 元 响应 不 一 致 性 等 均 影 响 测试 结果 可 靠 性 。 


所 以 除了 目标 图 像 Cimg 


) 外 ， 另 需 拍 摄 本 底 Chias) MFH lat) 图 像 进行 本 


底 消 除 和 平 场 改 正 来 消除 以 上 影响 ， 处 理 算 法 公式 为 


img — bias 
flat —bias 


基于 该 算法 ， 数 据 预 处 到 


imgr = 


C1) 


流程 简 述 为 


(1) 合并 多 张 bias 图 像 获得 本 底 文件 biasm. 
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(2) flat 文件 减 去 biasm 文件 ， 获 得 flatb 文件 。 


(3) 合并 多 张 flatb 文件 ， 获 得 平 场 文件 flatm。 


(4)img 文件 减 biasm 文件 , 再 除 以 flatm 文件 , 获得 真实 流量 数据 图 像 imgr。 


1.1.2 孔径 测 光 


对 于 图 像 中 较 小 的 天 体 ， 通 常 采 用 孔径 测 光 法 确定 其 亮度 。 算 法 包括 3 步 : 


d) 背景 天 光 确 定 和 消去 ; (2) 目标 检测 ; (3) 孔径 测 光 。 我 们 使 用 文 [11] 提 


出 的 算法 SEP (Source Extraction and Photometry) 完成 孔径 测 光 。 
由 于 大 气 存在 散射 现象 ， 散 射 光 入 入 CCD ， 成 为 数据 imgr 中 的 天 光 背 景 。 
SEP 将 图 像 分 割 为 矩形 块 ， 并 计算 每 个 矩形 块 中 所 有 像素 的 中 值 和 标准 差 等 。 对 


其 中 像素 点 应 用 中 值 滤波 ， 删 除 偏离 中 值 3 倍 标准 差 的 像素 点 ， 重 复 以 上 步 又 


直到 无 像素 点 可 删除 , 最 终 获 得 天 光 数 据 。 然 后 在 imgr 中 减 去 这 部 分 天 光 数 据 。 


SEP 使 用 一 种 多 阔 值 检测 算法 解决 目标 检测 问题 。 首先 设置 最 低 提 取 阔 值 找 


出 候选 点 , 按 提取 闵 值 和 峰值 划分 30 “SBOE, 基于 该 闵 值 把 图 像 构造 成 树 形 图 


= 


然后 从 树枝 开始 ,向 树干 进行 遍历 分 析 相 邻 点 是 否 属于 相同 的 目标 星 ， 从 而 


目标 星 的 位 置 和 形状 。 


基于 已 得 目标 位 置 , SEP 使 用 二 维 高 斯 曲面 拟 合 能 包括 目标 星 所 有 像素 点 的 


椭圆 孔径 ， 将 孔径 内 像素 点 亮度 积分 ， 获 得 imgr 文件 中 星 的 亮度 、 位 置 和 形状 


等 信息 ， 用 于 下 一 步 的 数据 提取 。 


1.2 目标 匹配 与 数据 提取 


由 于 望远镜 的 跟踪 误差 , 不 同时 刻 拍 摄 的 视 场 会 有 侦 移 ， 使 得 目标 星 在 图 像 


中 不 断 平移 。 目前 有 多 种 算法 用 于 匹配 星 图 目标 星 位 置 ， 如 基于 三 角形 "7 外、 基 


于 径 向 和 环 向 特征 "和 基于 向 量 特征 等 算法 。 我 们 对 现 有 的 算法 进行 比较 ， 
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结果 表明 ， 基 于 向 量 特征 的 算法 效率 较 高 。 但 星 图 偏 移 较 大 的 情况 下 ,存在 匹配 


失败 的 可 能 。 此 外 小 行星 会 在 星 图 中 不 断 移动 , 现 有 的 算法 无 法 自动 跟踪 小 行星 。 


本 文 提出 一 种 可 自动 追踪 目标 的 基于 向 量 特征 星 图 匹配 算法 。 
(1) 选择 一 张 星 图 作为 模板 图 像 ， 首 先 在 图 中 指定 待 提取 数据 的 目标 星 ， 


记录 目标 星 的 原始 坐标 值 , 然后 基于 该 模板 图 像 生 成 定位 向 量 特征 集合 。 生 成 方 


法 如 下 : 将 模板 图 像 中 心 区 域 星 排除 己 指 定 的 目标 星 外 的 其 余 星 按 亮 度 排序 ， 选 


前 mH 个 星 构成 定位 星 集合 ， 以 最 接近 该 集合 中 心 的 星 So 作为 定位 参考 星 ， 计 


算 其 余 m TEN $0 的 向 量 ， 构 成 定位 向 量 特征 集合 R= (5.7.07) « 


DEH R 匹配 其 余 星 图 文件 , 方法 如 下 : 将 待 匹 配 图 像 中 星 按 亮 度 排序 ， 
并 取 前 n 颗 得 待 匹 配 星 集 合 S ={51,5.…s,}。 设 si 为 参考 星 % ， 根 据 尺 计算 出 其 


R m 颗 定 位 星 的 坐标 ， 从 待 匹 配 图 像 中 寻找 m UE ce ETE, PRE ALE NT 


坐标 在 图 像 外 ， 则 假设 该 星 存在 ， 继 续 匹 配 其 余 星 


(3) WR m 颗 定 位 星 均 存在 ,说 明 匹 配 成 功 ，5, 相 对 于 So 的 偏 移 量 即 为 待 


匹配 星 图 文件 相对 于 模板 图 像 的 偏 移 量 。 根 据 该 偏 移 量 得 待 还 配 文件 中 目标 星 的 


= 新 位 置 ， 提 取 目 标 星 的 亮度 信息 。 


(4) 根据 目标 星 的 新 位 置 和 偏 移 量 ， 更 新 目标 星 的 原始 坐标 值 。 


(5) 如 未 能 发 现 全 部 m 颗 星 ， 说 明 匹 配 不 成 功 ，s1 并 非 参 考 星 5$,， 尝 试 假 


设 ss 星 为 参考 星 ， 依 次 类 推 。 如 果 n 颗 星 均 未 匹配 成 功 则 该 文件 匹配 失败 。 


相 比 于 传统 的 基于 向 量 特征 的 星 图 匹配 算法 , 本 文 算法 能 够 自动 处 理 图 像 偏 


移 太 大 导致 部 分 比较 星 偏 移 出 图 像 的 问题 


2 数据 实验 


天 文 数据 通常 使 用 FITS (Flexible image Transport System) 格 式 保 存 ， 需 要 用 对 应 的 程 
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序 包 分 析 和 读 取 ， 本 文 分别 以 双星 MMCom 和 小 行星 52Eur 为 例 ， 使 用 Python 


语言 中 astropy，numpy，sep，pandas 和 cv2 包 完 成 全 部 数据 处 理 流程 。 


(a) 原 始 图 像 (b) 预 处 理 后 图 像 
(a)before preprocessing (b)after preprocessing 


图 1 预 处 理 前 后 的 MMCom 数据 图 像 
Fig 1. The image of MMCom before and after preprocessing 
a) 观测 数据 目录 通常 包括 bias*.fits 〈 本 底 文件 )、flat*.fits 〈 平 场 文件 ) 


和 obj*.fits (目标 文件 ) 3 种 ， 我 们 对 文件 名 进行 扫描 ， 分 析 和 记录 下 完整 文件 
名 ， 然 后 自动 完成 数据 预 处 理工 作 。 (abefore preprocessing 
(b)after preprocessing 

1 (2) 和 (b) 分 别 为 MMCom 的 原始 数据 图 像 和 经 过 预 处 理 后 的 数据 图 像 。 


由 图 1 可 见 ， 经 过 预 处 理 后 ， 图 像 亮度 不 均匀 现象 消失 。 


(a) MMCom (b)52Eur 
图 2 全 部 测 光 结果 


Fig 2. photometry results 


《2) 使 用 SEP 对 预 处 理 后 的 图 像 进 行 孔径 测 光 ， 获 得 目标 位 置 与 大 小 信息 


如 图 2， 亮 度 、 位 置 和 孔径 等 信息 存 入 单独 的 CSV 文件 中 ， 供 后 续 流 程 使 用 。 


(a) MMCom (b)52Eur 


图 3 目标 星 及 向 量 特征 


Fig 3. object star and vector features 


(3) 取 一 张 星 图 作为 模板 图 ， 指 定 需 提取 的 目标 位 置 。 如 图 3， 其 中 红色 


圈 是 手工 指定 的 两 个 目标 ，1 号 为 目标 星 ，2 号 为 比较 星 。 图 3(b) 的 52Eur 不 断 


在 星 图 中 移动 。 随 后 程序 自动 选取 目标 以 外 最 亮 的 5 颗 星 用 于 定位 ， 其 中 RR 为 


Chinay in Ade HAI 
ChinaxIV 合 人 期刊 


参考 星 ，4 条 黄 线 是 用 于 星 图 匹配 的 向 量 集合 。 根 据 该 集合 ， 程 序 在 其 余 星 图 中 
自动 匹配 和 奶 踪 目标 并 提取 数据 。 


(4) 为 了 验证 所 得 数据 ， 我 们 用 IRAP 对 数据 进行 测 光 ， 并 将 两 种 方法 所 


得 结果 绘制 如 (a)MMCom (b)52Eur 


4。 从 图 中 可 以 看 出 ， 两 种 方法 测试 出 来 的 结果 基本 一 致 ， 其 中 MMCom 


和 52Eur 测 光 结果 最 大 偏差 分 别 为 0.04mag 和 0.043mag， 标 准 偏差 分 别 为 土 


0.005mag Fl +0.007mag. 
0.1 0.20 
024 N fa 
034 ff $ 
3 a r- & 
z 04 2 0.10 入 
:| r ` ] 3 So i F 出 
0.054 k z i 
074 r r ba 0.00 - 
9606.2 9606.3 9606.4 7419.1 7419.2 7419.3 
HJD(2450000.0+) HID(2450000.0 |) 
(a) MMCom (b)52Eur 


图 4 自动 测 光 与 IRAF 测 光 结果 比较 


Fig 4. comparison of photometry results by the program in this paper and by IRAF 


3 结论 与 展望 


(S) IRAF 虽然 功能 强大 ， 但 使 用 和 安装 流程 复杂 ， 且 二 次 开发 难度 较 大 。 本 文 


采用 Python 编写 一 套 全 自动 测 光 程序 ， 适 用 于 Windows, Linux 和 Mac 等 操作 


系统 ， 完 成 从 目录 扫描 分 析 文 件 类 型 ， 到 提取 最 终 数据 全 部 功能 ， 测 试 结果 与 


IRAF 一 致 , 具有 运行 效率 高 、 简 单 易 用 等 特点 ,希望 能 与 国内 天 文 观测 者 交流 ， 


进一步 完善 程序 功能 。 


目前 ， 淮 北 师范 大 学 天 文 研究 团队 于 2015 年 购置 30cm 望远镜 ， 最 佳 观 测 
星 等 约 为 7-10mag， 发 现 一 些 大 的 耀 发 说 ，2022 年 又 购置 50cm 平面 波 望远镜 ， 


处 于 调试 和 试 运行 阶段 。 我 们 计划 对 这 两 全 望远镜 (放置 地 点 为 国家 天 文 台 兴隆 
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基地 GWAC 组 ) 进行 远程 观测 改造 ， 将 获得 海量 的 测 光 数据 。 预 研一 套 自 动 测 


光 程 序 ， 这 为 测 光 数 据 的 实时 处 理 提供 了 可 靠 的 技术 保障 ; 同时 地 方 高 校 的 天 文 


发 展 提供 了 数据 平台 ,也 为 学 校 相 关 专 业 如 大 数据 、 光 电信 息 技术 专业 的 发 展 提 
供 支 持 。 
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Design and Implementation of Automatic Photometry Program Based on Python 
Yuan Huiyu', Zhao Juan’, Yang Yuangui’* 

(1.College of Information, Huaibei Normal University, 235000, Huaibei, Anhui, 
China 

2. Huaibei Normal University, 235000, Huaibei, Anhui, China) 
Abstract: This paper designs an automatic photometry program based on Python. 
Firstly, it scans the file name to achieve file classification; Then, it preforms data 
inspection, preprocessing, aperture photometry, and save the photometric results; 
Finally, after the observer specifying the target star and comparison star, the program 
extracts the corresponding magnitude and observation time from the photometric 
results to obtain the light curve based on the differential photometric method. In order 
to solve the problem of image matching in data extraction, an improved image 
matching algorithm based on vector features is proposed to automatically track 
moving objects, which can match the image with large offset and track the moving 
asteroids. The program is used to reduce the observational data of binary and asteroid 
respectively. Comparing the photometric results of this program with those of IRAF, 


the maximum deviations of photometry results are 0.04 and 0.043 mag respectively, 
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and the standard deviations are +0.005 and 0.007 mag respectively. 
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